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ABSTRACT: 

The invention concerns a method of determining the position or change in position or state of 
motion of a receiver or receivers, the position or change in position or state of motion of which is 
or are not already known, in a network of transmission sources some or all of whose positions 
are known. At a first time, the relative offsets in time, phase, frequency, or derivatives thereof with 
respect to each other, or with respect to a reference source, of the signals received by a first 
receiver from a plurality of the transmission sources are measured. At a selected second time, 
the respective offsets of the signals received by the same or a second receiver, whose position or 
state of motion is not already known, from the plurality of the transmission sources are measured. 
Then the position of the receiver or receivers at the first or second times or the change in position 
of the first receiver between the first and second times or the state of motion of either receiver is 
calculated from the relationship between the first and second sets of relative offsets 
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(57) Abstract: The invention concerns 
a method of deteimining the position or 
change in position or state of motion of 
a receiver or receivers, the position or 
change in position or state of motion of 
which is or are not already known, in a 
network of transmission sources some 
or all of whose positions are known. 
At a first time, the relative offsets in 
time, phase, frequency, or derivatives 
thereof with respect to each other, or 
with respect to a reference source, of 
the signals received by a first receiver 
from a plurality of the transmission 
sources are measured. At a selected 
second time, the respective offsets of 
the signals received by the same or a 
second receiver, whose position or state 
of motion is not already known, from 
the plurality of the transmission sources 
are measured. Then the position of 
the receiver or receivers at the first or 
second times or the change in position 
of the first receiver between the first 
and second times or the state of motion 
of either receiver is calculated from 
the relationship between the first and 
second sets of relative offsets. 
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IMPROVEMENTS IN RADIO POSITIONING SYSTEMS 

The present invention relates generally to improvements in radio positioning systems 
and their methods of operation and in particular to use in non-synchronised transmitter 
5 networks without the need for additional monitoring receivers. 

EP-A-0 303 371, the contents of which are hereby incorporated by reference, describes 
a radio navigation and tracking system which makes use of independent radio 
transmitters set up for other purposes. The signals from each transmitter, taken 

10 individually, are received by two receiving stations, one at a fixed and known location, 
and the other mounted on the mobile object whose position is to be determined. A 
representation of the signals received at one receiving station is sent via a link to a 
processor at the other receiving station, where the received signals are compared to 
find their phase differences or time delays. Three such measurements, made on three 

15 widely spaced independent transmitters, are sufficient to determine the position of the 
mobile receiver in two dimensions, i.e. its position on the ground. The phase or time 
offset between the master oscillators in the two receivers is also determined. 

"CURSOR", as the system described in EP-A-0 303 371 is known, is a radio 
20 positioning system which can use the signals radiated by existing non-synchronised 
radio transmitters to locate the position of a portable receiver. Unlike some other 
systems which use the temporal coherence properties of networks of purpose-built 
synchronised transmitters, CURSOR makes use of the spatial coherence of the signals 
transmitted by single transmitters. In a further development (see EP-A-0 880 712 & 
25 WO-A-99/21028), the technology has been applied to find the position of a mobile 
phone handset in a GSM or other digital telephone system, and these are examples of 
an 'Enhanced Observed Time Difference' (E-OTD) method using the down-link 
signals radiated by the network of Base Transceiver Stations (BTS) of the telephone 
system. 

30 

In the digital mobile telephone application described in EP-A-0 880 712, the contents 
of which are hereby incorporated by reference, the signals from each BTS within 
range of the handset are received both by the handset itself and by a fixed nearby 
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receiver, the Location Measurement Unit (LMU), whose position is accurately known. 
Representations of the received signals are passed to a Mobile Location Centre (MLC) 
where they are compared in order to find the time difference between them. In Figure 
1 we show the geometry of a standard two-dimensional system. The origin of 
5 Cartesian co-ordinates x and y is centred on the LMU positioned at O. The orientation 
of the axes is immaterial, but may conveniently be set so that the y axis lies along the 
north-south local map grid. The handset, R, is at vector position r with respect to the 
LMU position O. A BTS, A, is shown at vector position a. 

10 Consider first the signals from BTS A. The time difference, A/ a , measured between 
the signals received at R and O is given by 

Af a =(|r-a|-|a|)/t> + £, 

15 where v is the speed of the radio waves, e is the clock time offset between the clocks 
in the receivers at R and O, and the vertical bars each side of a vector quantity denote 
that it is the magnitude of the vector which is used in the equation. The value of £ 
represents the synchronisation error between the measurements made by the two 
receivers. Similarly, may written for two other BTSs (B and C) at vector positions b 

20 and c (not shown): 

Af b =(|r-b|-|b|)/v+£, 

Af c =(|r-c|-|c|)/u + £. (1) 

25 

The values of Ar a , Af b , Af c , are measured by the methods disclosed in EP-A-0 880 712 
and the values of a, b, c, and v are known. Hence the equations (1) can be solved to 
find the position of the handset, r, together with the value of e . 

30 In WO-A-99/21028, the contents of which are hereby incorporated by reference, it is 
described how these same time offsets can be measured using locally-created 
templates in a GSM telephone system as follows. Suppose that the handset has 
recorded a short burst of the GSM signals from BTS A. Contained within that 

2 
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recording is the framing structure, synchronisation bursts and other * given' data (or 
predetermined values) which are a constant feature of those transmissions. The 
processor within the handset can create a matching template, based on the known 
structure of the network signals. Received signals can then be matched by the locally- 
5 generated template. When the template finds a match, the correlation peak at the 
position of best match corresponds to the time offset between the received signals and 
the local clock inside the handset. For the signals radiated by BTS A this measured 
time offset, Af a i, is given by 

10 Afai =(|r - a|)/t> + e\ , 

where Cfe is the time offset of the BTS transmissions and £\ is the time offset of the 
handset's internal clock, both relative to an imaginary universal Absolute 1 clock. The 
signals from BTSs B and C may also be measured in the same way, giving 



15 



A/bi =(|r - b|)/i> + CM- e u 



and 



20 



Ar c i =(|r - c|)/u + ot+ £\ . 



(2) 



The same measurements can also be made by the LMU, giving 



Ara2 =(|a|)/u + Oa+ €z , 



25 



Af b 2 =(|b|)/i> + £ 2 , 



and 



30 



A/ c2 =(|c|)/v + «c+ € 2 , 



(3) 



where e 2 is the time offset of the LMLTs internal clock relative to the same imaginary 
universal absolute clock. Subtracting equations 3 from equations 2 gives 
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Af a = Ar a i - A/a2= (|r - a| - |a|)/u + e , 



A/ b = Af b i - Af b 2= (|r - b| - \b\)lv + e t 



5 



and 



Af c = A/ci - Afc2= (|r - c| - |c|)/t> + e , 



(4) 



10 where e = £| - e 2 . It will be noted that equations 4 are just like equations 1, and can be 
solved in the same way to find the position of the handset, r, and the value of e . 

The methods described above measure time offsets. However, it is sometimes useful 
to measure phase offsets, frequency offsets, or derivatives of frequency offsets as 
15 described later. 

It will be apparent that the CURSOR method, in common with all other methods 
which use the signals from non-synchronised transmitters, requires a network of 
LMUs to be set up within the coverage area of the telephone system. These units act 

20 as reference points at which the unsynchronised signals radiated by the BTSs are 
measured for comparison with the same signals received by a handset. In another 
patent application (our reference MJB06733WO) filed simultaneously herewith, we 
show how the entire network of BTS can be covered using just one "virtual LMU" 
which acts as the service node for all LMU data. The present invention shows how the 

25 CURSOR method (or other E-OTD method) can be applied without the need of a 
network of real LMUs. 

According to the first aspect of the invention there is provided a method of 
determining the position or change in position or state of motion of a receiver or 
30 receivers, the position or change in position or state of motion of which is or are not 
already known, in a network of transmission sources some or all of whose positions 
are known, the method comprising 
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at a first time, measuring the relative offsets in time, phase, frequency, 
or derivatives thereof with respect to each other, or with respect to a 
reference source, of the signals received by a first receiver from a 
plurality of the transmission sources; 

at a selected second time, measuring the respective offsets of the 
signals received by the same or a second receiver whose position or 
state of motion is not already known from the plurality of the 
transmission sources; and 

calculating the position of the receiver or receivers at the first or second 
times or the change in position of the first receiver between the first and 
second times or the state of motion of either receiver from the 
relationship between the first and second sets of relative offsets. 

By "state of motion" is meant a function of the movement of the receiver, including its 
15 velocity, acceleration, rate of change of acceleration etc. 

The relative offsets in time, phase, frequency, or derivatives thereof, with respect to 
each other or with respect to a reference source, of the signals received by the first or 
second receivers from a plurality of the transmission sources may be represented by 
20 corresponding offsets or differences in the distances between the transmission sources 
and the first or second receivers. 

It should be understood that, if the second time is measured at a second receiver, the 
first and second times could be the same. 

25 

At selected additional times, the respective relative offsets in time, phase, frequency, 
or derivatives thereof with respect to each other or with respect to a reference source 
of the signals received by the first receiver, or by the second receiver, or another 
receiver or receivers whose positions, changes in position, or states of motion is or are 
30 not already known, from a plurality of the transmission sources may be measured; and 
the position, change in position, or state of motion thereof of any or all of the receivers 
may be determined using the measurements and known positions of the transmission 
sources. 



(a) 



5 (b) 



(c) 



10 
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In the above aspect, it should be understood that any of the selected times could be the 
same as any other of the selected times. 

5 The invention also includes a system for carrying out this method, which system may 
include receivers controlled by computers or microprocessors suitably programmed. 

Examples of the method and system according to the present invention will now be 
described with reference to the accompanying drawings, in which: 

10 

Figure 1 shows the geometry of a standard two-dimensional system; 

Figure 2 is a diagram representing the plane of the Earth's surface (assumed to be flat) 
near to a pair of BTS units; 

15 

Figure 3 illustrates further the geometry of a standard two-dimensional system; 

Figure 4 shows a digital mobile telephone network with two handsets and three BTSs 
incorporating a system according to the invention; 

20 

Figure 5 shows a digital mobile phone network with three handsets and five BTSs 
incorporating a system according to the invention; and 

Figure 6 shows a hypothetical GSM network of BTSs and handsets. 

25 

To illustrate one example, in which the transmitters and receivers are part of a digital 
mobile telephone network, suppose that the handset is at vector position r(*i) at time 
U. Equations (2) are then 



Afai('i) = (|r(fi) - a\)/v + Oa(r,)+ , 



Ar bl (r,) = (|r(r,) - b|)/t> + Ob(fi)+ 



6 
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Afci(fi) = (|r(r0 - c\)/v + Ofc(r0+ 



(5) 



where Oa(fi) is the time offset of the BTS transmissions from A, OfeC'i) the offset from 
B, Oc(t\) the offset from C, and £\(t\) is the time offset of the handsets internal clock, 
5 all measured at time t\ relative to an imaginary universal 'absolute* clock. Now, if the 
handset is at a different vector position r(f 2 ) at a later time r 2 , equations (2) are then 



Subtracting equations (6) from equations (5) gives 

15 

Afai(*i) - A/ a ,(r 2 ) = (|r(fi) - a| - |r(r 2 ) - a|)/o + (Oa(f,) - Oa(r 2 ))+ (Ci('i) - e } (t 2 )) , 

Ar b ,(r,) - Ar bl (/ 2 ) = (|r(;,) - b| - |r(* 2 ) - b\)/v + (atfr i) - ab(f 2 ))+ (£i(f 0 - £i(r 2 )), 

20 Afcfr) - A/ Cl (r 2 ) = (|r(r,) - c| - |r(r 2 ) - c|)/v + (ot(r,) - Ofc(f 2 ))+ - £,(f 2 )). (7) 

As has been noted above, the BTS transmissions in a GSM or other digital network are 
not usually synchronised. However, the BTSs are usually equipped with high-stability 
oscillators and the signals are often locked to a common reference timing source. 
25 Where this is the case, the signals keep a constant time offset with respect to each 
other, so £a('i) = £*(ti), £b('i) = £b(fe)« and = edti). (Even where this is not the 
case, the BTS signals are often sufficiently stable to allow an estimate to be made of 
these differences in a service node or other network element.) Writing Ar a = Ar a i(fi) - 
A/ a i(f 2 ), Ar b = Ar b i(/j) - Ar b i(f 2 ), A/ c = Ar c i(f,) - Af cl (r 2 ), and e = £,(r,) - £,(f 2 ), gives 



A/ a ,(f 2 ) = (|r</ 2 ) - a|)/t> + Oa(f 2 )+ £,(f 2 ) , 



10 



Afbifo) = (|r(r 2 ) - h\)lv + ab(/ 2 )+ £\(t 2 ). 



Ar cl (r 2 ) = (|r(r 2 ) - c|)/i> + ofc(/ 2 )+ 6,(r 2 ). 



(6) 



30 



A/ a = (|r(r,) - a | - |r(r 2 ) - a|)/u + e, 



A^^draO-bl-lrC^-bp/v+e, 

7 
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At c = (|r(*,) - c| - |r(r 2 ) - c\)/v +e. (8) 

Equations (8) can be understood graphically by reference to Figure 2. The diagram 
5 represents the plane of the Earth's surface (assumed to be flat) near to BTS units A 
and B. In the first of the equations, the term |r(fj) - a| represents the distance of the 
handset from A at time t u and the term |r(r 2 ) - a| represents its distance at time t 2 . The 
first of equations (8) can be rewritten as 

10 |r(f,) - a| - |r(* 2 ) - a) = vAt a -ve 9 

which therefore represents the loci of two concentric circles centred on A whose radii 
differ by vAt a - ve. These are marked 1 and 2 respectively in Figure 2. The above 
equation does not define the radius of either circle, but only the difference between 

15 them. By itself, this equation therefore does nothing to locate either of the points r(/j) 
or r(f 2 ). Also marked in Figure 2 are a second pair of concentric circles, 3 and 4, 
representing the second of the equations (8) and centred on B. Again, their radii are 
undefined by the equation, but the difference between their radii must be vAtb - V£. 
Point r(fi) must lie at one of the intersections of circles 1 and 3, and point r(r 2 ) must 

20 lie at one of the intersections of circles 2 and 4. Suppose that the value of £ is zero, i.e. 
that the handset's internal clock has kept perfect time between the measurements; then 
if we know the position of point r(fj), say P in Figure 2, we can deduce that point r(f 2 ) 
must be at Q, because the four circles are now fixed in space, and hence we have 
measured the position of the handset at the later time. In practice, we can't assume 

25 that £ is zero, so we must use all three of equations (8) to find Q given P. 

It will be noted there is also a second point of intersection, Q 1 , of the circles 2 and 4, 
and hence an ambiguity in the determination of r(f 2 ). If this ambiguity cannot be 
resolved by other means (for example, by knowing that the handset is being used by a 
30 pedestrian who could not have moved from P to Q' in time t 2 - t\) then four 
measurements involving four BTSs may be used. 



8 



WO 00/73814 



PCT/GBOO/02105 



One of the advantages of the present invention is that if a single handset makes 
measurements of the signals from a minimum of three geographically separate BTSs at 
two different times, then it is possible to determine the change in the position of the 
handset between these two times (without reference to a known starting point), 
5 provided that at least three of the BTSs are common to the two sets of measurements. 
If the position of the handset has been determined at some point previously, use of this 
technique can enable applications requiring autonomous navigation within the handset, 
i.e. navigation without further reference to the CURSOR or any other E-OTD system. 

10 An extension of the invention, described later, shows how the times t\ and tz can be 
brought arbitrarily close together so that the measurements produce estimates of rates 
of change, or derivatives, of the measured quantities. 

The above discussion shows how the position of a handset at a later time can be 
15 obtained from three CURSOR-like measurements on three geographically-dispersed 
BTSs made both at the later time and at an earlier time provided (a) that the position of 
the handset is known at the earlier time, and (b) that there is no relative drift between 
the signals radiated by the BTSs in the interval (or that such drift is known). The 
accuracy of the position determination depends both on the precision with which the 
20 first position is known and on the distance moved between the measurements. 
Equations (8) actually contain five unknown quantities: the x and y co-ordinates of 
each of the two points P and Q, and the unknown handset clock drift e. CURSOR-like 
measurements on each of five geographically-dispersed BTSs at both a first and a 
second position are therefore sufficient to determine both P and Q uniquely if the 
25 handset has moved appreciably between the two sets of measurements. 

Thus, the invention provides an E-OTD method of finding the position of a moving 
handset in an unsynchronised digital mobile phone network without the need of any 
LMUs at all. It may be especially useful for tracking the handset in a network which 
30 does not already have a location system. 

A further advantage of the present invention is the ability to use the known position of 
one handset to find the unknown position of a second handset. Suppose that each 

9 
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handset makes a measurement at about the same time using three geographically- 
dispersed BTSs A, B, C (see Figure 4). Equations 5 for the first handset, Hi, are then 

Ar ai (ri) = (|ri(*i) - a|)/v + Oa(fi)+ £i(f,) , 

5 

M,,(r,) = (|n(*,) - b\)/v + ob(ri)-f ei(fi), 

Afci('i) = (|ri(/,) - c|)/t> + ot(/,)+ e,(/i), (9) 

10 where ri(fj) is the (known) position of and the measurements are made at time t\. 
The same equations for the second handset, H 2 , are 

Afaite) = (|r 2 (r 2 ) - a|)/t> + OateH ^2) , 

15 A/ b ,(* 2 ) = (|r 2 (f 2 ) - b|)/v + ObteH ^('2), 

Afcifo) = (|r 2 (r 2 ) - c|)/v + ot(f 2 )+ eiihX (10) 

where r 2 (r 2 ) is the (unknown) position of H 2 and the measurements are made at time r 2 . 
20 Subtracting equations (10) from equations (9), noting that Oa(/i) = Oa(f 2 ), <*b(*i) = 
Obfo), and a^(ri) = a^(r 2 ), and writing A/ a = Ar a i(*i) - Ar a i(r 2 ), At b = Ar b i('i) - Ar b i(f 2 ), 
Ar c = Ar c i(fi) - Ar c i(f 2 ), and e = ,) - e 2 (r 2 ), we get 

Ar a = (|n(f 1) - a| - |r 2 (r 2 ) - *\)'v + e, 

25 

ATb = (W^)-N-W^- b IVi> + ft 

A/ c = (|ri(/i) - c| - |r 2 (/ 2 ) - c|)/u +£, (11) 

30 which are identical in form to equations (8). Hence, given ri(f a ), the position of Hj at 

time tu we can calculate r 2 (/ 2 ) which is the position of H 2 at time r 2 . (Again, 

measurements on four BTSs are needed to resolve the ambiguity in r 2 (r 2 ) if this cannot 

be resolved by other means.) In this case, t\ may be equal to / 2 . 

10 
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This same idea can be extended to many handsets. In equations (1 1), the vector ri(t\) 
refers to a first handset Hi whose position is known, and the vector r 2 (r 2 ) refers to any 
second handset H 2 - H n , perhaps representing any one of a large number. It is 
5 therefore possible to use, as a temporary measure, a first handset Hi at a known 
location as the "LMtT, offering the possibility of establishing an E-OTD service very 
quickly in a new area. 

This method can be used even in an area already covered by an LMU. Having used a 
10 normal E-OTD method to find the position of a handset H m the handset can then be 
tracked without further reference to the LMU. The calculation of the handset's 
position could be made at a service node (e.g. the mobile location centre; MLC), at 
another location, or in the handset itself, depending on the application. 

15 Another interesting application is to a system of many handsets Hi - H n in an area in 
which none of their positions is known. Provided that any pair sufficiently far apart 
can measure the same five BTSs A - E (see, by way of example, Figure 5), their 
positions can be determined using five equations like the three in (8). Otherwise, 
provided that separated pairs can measure at least three or four common BTSs, 

20 calculations based on many such measurements by pairs over a short period may be 
sufficient to determine the positions of all of them. 

Consider, for example, a set of three handsets, Hi, H 2 , and H 3 (see Figure 5). The 
relative positions of Hi and H 2 , H 2 and H 3 , and H 3 and Hi, may be computed using 
25 three or four BTS measurements since the constraints imposed by the equations on the 
positions of the vertices of the triangle HiH 2 H 3 may make the solution unique. In any 
case, well-known statistical techniques exist, such as those based on Bayes' theorem 
or a maximum likelihood method, which can be used to solve the equations. 

30 In the above discussions, we have made the assumption that the BTS network stability 
is good enough to be able to predict relative drifts between measurements. It is 
possible to apply the invention to a network in which the BTSs signals are free- 
running with respect to each other, i.e. one in which there is no synchronisation of any 



11 



WO 00/73814 



PCT/GB00/02105 



sort between the transmissions. In this case, at least two handsets are required which 
can receive the signals from the same five geometrically-dispersed BTSs at 
substantially the same time as each other. Equations (11) then become 

5 Ar a = (|n(fi) - a| - |r 2 (ri) - a|)/u + £, 

Ar b = (|n(f,) - b| - |r 2 (*i) - b|)/t> + £, 

Af c = (\ri(t x ) - c| - |r 2 (r!) - c\)/v + e, 

10 

A/^dn^o-dl-lnC/O-dD/o + e, 



Ar e = (Mr,) - e| - |r 2 (r,) - e\)/v +£, (12) 

15 where a, b, c, d, and e are the vector positions of the five BTSs, Ti(t\) is the vector 
position of the first handset H t and r 2 (fi) is the position of the second handset H„, both 
sets of measurements being made at the same time fi. This moment could be 
signalled, for example, by a particular element of the transmissions such as a particular 
frame number or a special signal. In equations (12), the transmission time offsets of 

20 the signals from the BTSs (the quantities Oa, 0Cb» otc etc.) cancel out and do not appear 
provided that the measurements are made by the two handsets sufficiently close 
together in time that the drifts are small enough to ignore. The five measurements in 
equations (12) are sufficient to locate the two handsets. Note that no LMU or 
equivalent is required in this case. 

25 

It is described in another copending PCT patent application (our ref: MJB06733WO), 
the contents of which are hereby incorporated by reference, how the measurements 
from a network of LMUs can be collected together into one list, the Virtual' LMU 
(VLMU) as if all the BTSs could have been measured by that single VLMU. An 
30 extension to this invention is to set up one or more VLMUs in a GSM or other 
transmission network which has no real LMUs at all. Timing measurements by 
handsets like those described above, especially in relation to equations (5) to (12), can 
be processed to provide a map of the receive time offsets of the signals from the BTSs, 
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which may be kept in the VLMUs and used in subsequent standard CURSOR or other 
E-OTD position calculations. One might also imagine a system in which a sparse 
coverage of real LMUs is augmented in one or more VLMUs by handset 
measurements to provide the same level of service as from a full network of real 
5 LMUs. 

The operation of an exemplary system incorporating the invention will now be 
described with the aid of simulated measurements. 

10 In Figure 6 is shown a hypothetical network of five BTSs B by the filled circles. Also 
shown are the positions of two handsets H by the open squares, a distance of about 2.2 
km apart. A computer program (which in practice could be running in the 
microprocessor of one or both of the handsets) was used to generate a list of the timing 
offset measurements which would have been made by each of the handsets H, 

15 including the addition of a random noise element equivalent to 50 metres rms to 
simulate real conditions more closely. Each BTS B was assigned an arbitrary and 
different signal timing offset, a, selected at random with respect to a universal 
•absolute 1 clock. Each handset H was also assigned a random value for the 
synchronisation error, £, of its internal clock. The distances from each handset to 

20 each BTS were then used to compile a list of the time offsets of the signals received by 
each handset relative to its internal clock. The two lists of timing offsets were then 
processed in a computer programmed to solve the five equations (12) above. The 
results are shown in the figure by the crosses X. The errors in eastings and northings 
for the handset towards the lower left-hand corner were -95 and -59 metres 

25 respectively, and for the handset near the upper right-hand corner the errors were -100 
and +30 metres. The ellipses F surrounding the marked positions represent the regions 
of confidence within which the true position was predicted to lie at the 67 per cent 
level. The semi-major axes were of length 160 and 240 metres. 

30 It is interesting to note that, in general, there is a significant correlation between the 
errors in the positions generated by this method, with the result that the distance from 
one handset to another is known more accurately. In the example above, the error in 
the measured distance, actually 2.2 km, was less than 90 metres. 

13 
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A second example of a system according to the invention will now be described in 
which various rates of change, or derivatives, of the measured quantities are used. 

5 In the first of equations 7 the case was considered of a single handset making 
measurements of the signals from a BTS, A, of the observed time offset, A/ a i(*i), when 
at position r(fj) at time tu and of the observed time offset, A/ a i(/ 2 ), when at position 
r(/ 2 ) at time f 2 . This equation is reproduced here: 



10 ATai(fi) - A/ al (r 2 ) = (|r(r,) - a| - |r(r 2 ) - a|)/u + (Oa(/i) - a£t 2 ))+ (£,(*,) - f^fo)). 

In what follows, it is helpful to adopt the following changes to our nomenclature: 
/] becomes t 

t 2 becomes t + At 

15 A*ai(0 becomes T a (f) 

T(r + At) - T(r) becomes AT a (f) 
£ { becomes e . 

The equation then becomes 

20 

AT a (/) = (|r(H-Ar) - a| - |r(r) - a|)/u + a^r+Af) - Oa(f) + £(t+At) - £ (f). (13) 



The assumption is made once again that the stability of the BTS transmissions in the 
network is such that the values of a do not change significantly over the period of the 
25 measurements, At, (which can be made arbitrarily short). Thus Oa(H-Ar) = Oa(0* and 
equation 13 can be simplified. The complete set for five BTSs then becomes 



AT a (r) = (|r(/+Af) - a| - |r(r) - a|)/u + £(t+At) - £(t) 
30 AT b (/) = (|r(r+A0 - b| - |r(<) - b\)/v + e(r+A/) - £(t) 

AT c (r) = (|r(/+Ar) - c| - |r(r) - c\)/v + £(/+A/) - e(t) 



14 
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AT d (r) = (|r(f+Ar) - d| - |r(r) - d|)/t> + e(f+Ar) - e(t) 



AT e (0 = (|r(f+Ar) - e| - |r(/) - e\)/v + e(t+At) - e (t). 



(14) 



5 These equations can be solved in precisely the same manner as equations 12 to give 
both r(t+At) and r(f) without reference to an LMU network. As Ar is allowed to tend 
to zero, so rates of change, or derivatives, of the quantities are found as follows: 



10 



dt 



dt 



b _ 



dt 



dt 



d_ _ 



dt 



V COS0 a /v + — 
1 a dt 



y\cos6 b /v- 



d£ 
dt 



v|cos0 c /u + — 
dt 

vlcos0 d /t>+— 
1 dt 

v|cos0 e /u+— , 
dt 



(15) 



15 where v is the velocity of the handset (v = dr I dt), 6 is the angle between the velocity 
vector and the vector from the BTS to the handset (see Figure 3), and de I dt is the rate 
of change of the handset's internal clock, i.e. the frequency offset of the handset's 
crystal oscillator from its nominal value. These equations are mathematically the 
same as equations 14 and can be solved to give both the handset's position, r, and 

20 velocity, v, without reference to an LMU network. 

If it is also assumed that the handset is travelling at a constant velocity, then further 
derivatives can be evaluated as follows: 



25 



v „ a = , 1 ' , sin fl a , 



dt' 

d% 
dt' 



r-a 



= -3- 



r-a 



~cos0 a sin 2 0 a , 



(16) 
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in which the assumption has been made that the value of d 2 e I dt 2 is negligible. 

In many circumstances, measurements of only some of these parameters might be 
available. Nevertheless, they can be included with other E-OTD calculations to 
5 improve the estimates of the handset's position. However, if measurements are made 
of all of these parameters then both the distance between the BTS and the handset, and 
the handset's velocity, can be determined as functions of del dt as follows 



10 



1 [ dt dt 1 1 d? ) y dt 2 1 d? > 



de 
dt 



II a dT * de 



Jfc_dT^ 
dt dt 



(17) 



dt 3 



15 



20 



If measurements are made on three or more geographically disparate BTSs, then one 
can solve directly for r, de I dt and v without any reference to an LMU network. 

There are circumstances, however, when the equations cannot be evaluated. These 
occur when d^T* / dT 3 is zero. Even in these cases information about r and v can still 
be obtained as there are only four possible reasons for this zero value (in the absence 
of noise): 

1. If the handset is stationary, ie. v = 0, then both d 2 T a I dt 2 and c?T a I dt 3 are 



zero, so 



dT a _ de 



dt dt 

2. If 0 a = ±90° then tPT* I dt 3 will be zero and hence 
d 2 T. 



dt' 



dT a _ de 
|r-a| v' and dt dt 
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3. If the handset is travelling directly towards the BTS, i.e. 0 a = 0°, then d 2 T a I dl 2 

and <?T a / dt* will both be zero, whilst = M 

dt v dt 

4. If the handset is travelling directly away from the BTS, i.e. 0 a = 180°, then both 

d 2 T a I dT 2 and d*T a I dt 3 will be zero, whilst ^ = ztl + — . 

dt v dt 



If the handset passes directly underneath the BTS while making measurements, i.e. 

r = a, then there is a discontinuity in f^ZL. If such a discontinuity is observed by a 

dt 

handset then its position is immediately known through a measurement from just the 
single BTS. 



10 



When handsets are used in a digital telephone network, such as a GSM system, the 
handset timing is made to advance with respect to the signals received from the 
'serving BTS* (ie the BTS with which the handset is in communication for telephonic 
purposes), so as to cause the handset signals to be received by the serving BTS in 
15 synchronism with its own signals. This is often known as 'timing advance'. In 
practice, this means that the handset synchronises its own internal clocks to the signals 
received from the serving BTS, say BTS A. This synchronisation alters equations 16 
to have the following forms: 



20 = Iv|(cos0 b - cos0 a ) 

dt 



dt 1 



v — ^ = |v| 



'sin 2 0 b sin 2 ©, > 



I k-b| 



i r - a i , 



dt 3 1 1 



cos 6 b sin 2 fl b cos fl a sin 2 0, ^ 



(18) 



Measurements of some or all of these values also can be combined with the other E- 
25 OTD measurements to give additional constraints and hence better accuracy of the 
predicted position and velocity. For example, if all of the above parameters are 
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measured for both BTSs B and C, then as before the equations can be solved without 
reference to an LMU network to give the position and velocity of the handset. 

In the foregoing example of a system of the invention, measurements of timing offsets 
5 have been used, for example to estimate the rate of change of a timing offset. This is 
equivalent to measuring the frequency offset of the signals from that BTS, and this 
may be achieved directly by measuring and analysing the base-band signals from the 
receiver. For example, the digitised I and Q data samples corresponding to, say, a 
synchronisation burst in a GSM system (whose form is known in advance) could be 
10 transformed into the frequency domain using a well-known technique such as a fast 
Fourier transform, and the frequency of identified features compared with expected 
values to obtain the frequency offset of the signals. 

This technique has been simulated in a computer program to test the effect of the 
15 signal to noise ratio on the accuracy with which the speed of the handset can be 
measured by this method. The accuracy depends also on the number of samples, N, in 
the template used to match the fequency spectrum. The results are as follows: 



S/N ratio 


N = 64 


N = 128 


N = 256 










10 


38.4 


12.8 


4.8 


30 


12.8 


4.3 


1.6 


100 


3.8 


1.3 


0.5 











20 This table gives the error in kilometres per hour associated with corresponding values 
of S/N ratio and template length. 
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CLAIMS 

1 . A method of determining the position or change in position or state of motion 
of a receiver or receivers, the position or change in position or state of motion of 
5 which is or are not already known, in a network of transmission sources some or all of 
whose positions are known, the method comprising 

(a) at a first time, measuring the relative offsets in time, phase, frequency, or 
derivatives thereof with respect to each other, or with respect to a reference 
source, of the signals received by a first receiver from a plurality of the 

10 transmission sources; 

(b) at a selected second time, measuring the respective offsets of the signals 
received by the same or a second receiver, whose position or state of 
motion is not already known, from the plurality of the transmission 
sources; and 

15 (c) calculating the position of the receiver or receivers at the first or second 

times or the change in position of the first receiver between the first and 
second times or the state of motion of either receiver from the relationship 
between the first and second sets of relative offsets. 

20 2. A method according to claim 1, wherein the second time is measured at a 
second receiver and the first and second times are the same. 

3. A method according to claim 1 or claim 2, wherein at selected additional times, 
the respective relative offsets in time, phase, frequency, or derivatives thereof with 

25 respect to each other or with respect to a reference source of the signals received by 
the first receiver, or by the second receiver, or another receiver or receivers whose 
position, change in position, or state of motion is or are not already known, from a 
plurality of the transmission sources is measured; and the position, change in position, 
or state of motion thereof of any or all of the receivers is determined using said 

30 measurements and known positions of the transmission sources. 

4. A method according to claim 3, wherein any of the selected times is the same 
as any other of the selected times. 
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5. A method according to any of claims 1 to 4, wherein said receivers are 
handsets of a GSM system and wherein the or each handset has its timing advanced to 
be in synchronism with a serving BTS. 

6. A method according to claim 5, which further includes 

(d) measuring, at a first time, the relative offsets in time, phase, frequency, or 
derivatives thereof, with respect to each other or with respect to a reference 
source of the signals received by a first handset from a non -serving BTS; 

(e) measuring, at a selected second later time, the relative offsets in time, 
phase, frequency, or derivatives thereof respectively with respect to each 
other or with respect to a reference source of the signals received by the 
same or a second handset from the non-serving BTS; 

(f) calculating the position of either handset or the change in position of the 
first handset between the first and second times or the state of motion of 
eith er handset atjhe Jirst^or second times from the relationship between the 
first and second sets of data. 

7. A system for determining the position or change in position or state of motion 
of a receiver or receivers, the position or change in position or state of motion of 
which is or are not already known, in a network of transmission sources some or all of 
whose positions are known, the system comprising 

(a) means for measuring, at a first time, the relative offsets in time, phase, 
frequency, or derivatives thereof, with respect to each other or with respect 
to a reference source of the signals received by a first receiver from a 
plurality of the transmission sources; 

(b) means for measuring, at a selected second time, the respective offsets of the 
signals received by the same or a second receiver, whose position or state 
of motion is not already known, from the plurality of the transmission 
sources; 

(c) means for calculating the position of the receiver or receivers at the first or 
second times or the change in position of the first receiver between the first 
and second times or the state of motion of either receiver from the 
relationship between the first and second sets of relative offsets. 
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8. A system according to claim 8, further comprising: 

means for measuring, at selected additional times, the respective relative 
offsets in time, phase, frequency, or derivatives thereof with respect to each other or 
5 with respect to a reference source of the signals received by the first receiver, or by the 
second receiver, or another receiver or receivers whose position, change in position, or 
state of motion are not already known, from a plurality of the transmission sources; 
and 

means for calculating the position, change in position, or state of motion of any 
10 or all of the receivers using said measurements and known positions of the 
transmission sources. 

9. A system according to claim 7 or claim 8, wherein said measuring means and 
said calculating means include a microprocessor in one or more of said receivers. 

15 

10. A system according to any of claims 7 to 9, which comprises a component or 
components of a GSM system. 

11. A system according to claim 10, which includes one or more handsets, the or 
20 each of which has its timing advanced to be in synchronism with a serving BTS. 

12. A system according to claim 11, which further includes 

(d) means for measuring, at a first time, the relative offsets in time, phase, 
frequency, or derivatives thereof, with respect to each other or with respect 

25 to a reference source of the signals received by a first handset from a non- 

serving BTS; 

(e) means for measuring, at a selected second time, the relative offsets in time, 
phase, frequency, or derivatives thereof respectively with respect to each 
other or with respect to a reference source of the signals received by the 

30 same or a second handset from the non-serving BTS; 

(f) means for calculating the position of either handset or the change in 
position of the first handset between the first and second times or the state 
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of motion of either handset at the first or second times from the relationship 
between the first and second sets of data. 
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Figure 2 
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Figure 5 
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